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Q\ ', Abstract. The Lin and Pringle model (1987) of galactic disk formation 

^/^ • postulates that if star formation proceeds on the same timescale as the 
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Introduction 



viscous redistribution of mass and angular momentum in disk galaxies, 
^^ ' then the stars attain an exponential density profile. Their claim is that 

this result holds generally: regardless of the disk galaxy's initial gas and 
dark matter distribution and independent of the nature of the viscous 
processes acting in the disk. We present new results from a set of 2D 
^^ I hydro-simulations which investigate their analytic result. 

m ' 

o. 

^^ . A popular explanation for exponential stellar light profiles of both high and 

^^ ' low surface brightness disk galaxies is that stars form in situ in an already 

^^ . exponential gaseous disk (Mo, Mao & White 1998 and references therein). The 

Qh' origin of the exponential gaseous disk is in turn explained by the idea that all 

Q . gas destined to settle into the disk originated in a rigidly rotating spherical 

H I protocloud which conserved its specific angular momentum not only during the 

^ . process of protogalactic collapse and settling into centrifugal equilibrium but 

also throughout the star formation process (Mestel 1963; Gunn 1982). Given the 
hierarchical picture of galaxy formation which makes such idealized conditions 
unlikely, a more appealing idea is that stars acquire an exponential stellar profile, 
;2j ■ regardless of initial conditions prior to star formation, if star formation proceeds 

5^ I on a viscous timescale (Silk &: Norman 1981; Lin &: Pringle 1987). 

We extend analytical work done by previous authors by performing numer- 
ical hydrodynamical simulations with the BGK (Bhatnagar-Gross-Krook) code 
(Prendergast & Xu 1993; Slyz & Prendergast 1999). Because this code solves 
a model to the collisional Boltzmann equation, it incorporates dissipation natu- 
rally. Our simulations start with a gaseous disk of arbitrary initial density profile 
in inviscid centrifugal equilibrium in a static dark matter halo also of arbitrary 
profile. On a 201 by 201 evenly spaced Cartesian grid, we follow the evolution 
of disks with Milky Way type masses and sizes for about 12 Gyr. We assume 
an isothermal equation of state and we neglect self-gravity and stellar feedback 
in this contribution. 

2. Results 

For different initial gas and dark matter profiles we obtain disks with approxi- 
mately exponential stellar profiles as long as t^ ^ tvisc (Figs. 1 and 2). We find 

1 



Figure 1. Comparison of stel- 
lar profiles after ~ 12 Gyr for 
different values of /3 = tv,/tvisc- 
The gravitational potential, <1>, is 
ivo2ln(ro2 + r^) {vq = 220 km/s, 
vq = 2.56 kpc). The initial gas 
profile is nearly frozen into the 
stellar profile for small /3 whereas 
for large /? too much viscous evo- 
lution occurs and the resulting 
stellar profile is more of a power 
law than an exponential. 
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Figure 2. Comparison of stellar profiles after ~ 12 Gyr for (3 = 1 and 
for different rotation curves (rg in ^ is 1.3 kpc on the left and 10.2 kpc 
on the right). The initial gas profile is uniform with S = 39.8 Mq/pc^. 



that the stellar exponential scalelength is approximately half the initial half mass 
radius for the gas. Hence in so far as the disks all end up with exponential stel- 
lar profiles, they forget their initial conditions but the exponential scalelength 
remembers the degree of central concentration of the initial gas distribution. 
In addition to obtaining exponential stellar profiles, we also find gas and star 
fractions after ~ 12 Gyr that are consistent with observations without having to 
invoke an efficiency parameter for star formation. The exact link between star 
formation and viscosity remains a question. However our results suggest that 
the viscous timescale is indeed the natural timescale for star formation. 
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